clear;close all;
% file = '.\S110deg.txt';
% file = '.\S1148deg.txt';
% file = '.\S1190deg.txt';
Data = readCST(file, 4, 5);
Fre = Data(:, 1);
S11 = Data(:, 2) + 1j * Data(:, 3);
Z = sqrt((1 + S11).^2 ./ (1 - S11).^2);
SAVE = [real(Z),imag(Z)];
YW_Plot_Fast({Fre, real(Z), Fre, imag(Z)}, {"*-", "v:"}, ["Real", "Imag"], ["Frequency (THz)", "Impedance"]);
